home *** CD-ROM | disk | FTP | other *** search
/ Libris Britannia 4 / science library(b).zip / science library(b) / HAMRADIO / APRS400.ZIP / MAPFIX.BAS < prev    next >
BASIC Source File  |  1994-02-14  |  54KB  |  1,170 lines
  1. REM MAPFIX.bas PROGRAM.  SEE EXPLAINATION BELOW
  2. REM
  3. Ver$ = "4.00"
  4. MaxNumMAPS = 120'was 99 Current maximum number of maps loaded by APRS
  5. MaxNumPoints = 1500 'was 1000
  6. MaxNumLABELS = 99 'was 99
  7. MaxNumLines = 900
  8. REM $DYNAMIC
  9. GOTO BEGIN
  10.  
  11. Info: COLOR 15, 4: CLS
  12.    PRINT " MAPFIX.bas VERSION "; Ver$; " PROGRAM FOR FIXING APRS MAPS": PRINT
  13.    PRINT " This program has evolved significantly.  305 added a TRIM command to eliminate"
  14.    PRINT " all points outside an area.  3.07B adds W7KKE's GPS track history overlays and"
  15.    PRINT " digitizer code for using a drawing tablet.  308 fixed a filename bug.  Now in"
  16.    PRINT " 312, it will take USGS data from the 2,000,000:1 CD ROM and build APRS maps!"
  17.    PRINT " BUT NOT EASILY!  For a 50 mile map, there are 5 MBytes of USGS data which must"
  18.    PRINT " be filtered to get down to the 10K APRS size for the same area! The final steps"
  19.    PRINT " are all manual and take almost as long as developing an APRS map from scratch!"
  20.    PRINT
  21.    PRINT " CAUTION, THIS PROGRAM IS NOT PERFECT... KEEP BACKUPS!  Do a little at a time!"
  22.    PRINT
  23.    PRINT " Although MAPFIX has commands to make map modifications easier, a text EDITOR is"
  24.    PRINT " still useful for whole scale rearranging of points and features in a map file. "
  25.    PRINT
  26.    PRINT " MAPFIX uses two cursors.  The normal yellow APRS cursor, and a White MapPoint"
  27.    PRINT " which will be the next point to be processed.  ALT Keys allow you to MOVE the"
  28.    PRINT " MapPoint to the cursor, ADD a new point at the cursor, or DELETE the MapPoint."
  29.    PRINT " MAPFIX.bas shows you Deg/Min, Decimal, and APRS values of the cursor position."
  30.    PRINT
  31.    PRINT " ALSO NOTE THAT THE LIMITS IN APRS ARE 1000 POINTS, 99 FEATURES, and 99 LABELS!"
  32.    PRINT " If you need more points, features or labels, begin another map.";
  33.    LOCATE 25, 1: PRINT " HIT ANY KEY to continue...";
  34.    GOSUB GetChar: BEEP
  35.   
  36. Info2:   COLOR 15, 8: CLS
  37.    PRINT " PAGE 2 INSTRUCTIONS:  More about new features in version "; Ver$
  38.    PRINT
  39.    PRINT " With this new MAPFIX.bas, you can not only modify features by moving, adding"
  40.    PRINT " or deleting POINTS, but you can now add and kill FEATURES too, ie: roads,"
  41.    PRINT " rivers, borders, etc, from within the program.  In addition to the new KILL"
  42.    PRINT " and TRIM commands, the G key will move the cursor (GOTO) to the MapPoint,"
  43.    PRINT " and the F key will FIND the MapPoint nearest the cursor location.  If the"
  44.    PRINT " MapPointer and FeatureName get out of sequence, the RESET command may fix"
  45.    PRINT " them, but you should save the file immediately and check it with an editor."
  46.    PRINT
  47.    PRINT " I find the capability to delete points very useful when making larger area"
  48.    PRINT " maps from several smaller detail maps.  First, I run MAPCNVRT.bas to convert"
  49.    PRINT " all of the smaller maps to new temporary files with the new origin of the new"
  50.    PRINT " larger map.  Then I use the KILL command in MAPFIX to eliminate all minor "
  51.    PRINT " roads, features and labels and then the DELETE POINT command to remove all"
  52.    PRINT " inconsequential minor points from the roads that will not be needed"
  53.    PRINT " at the larger scale.  Then I use the editor to combine all of the points and"
  54.    PRINT " labels into the new file."
  55.    PRINT
  56.    PRINT " A new MAPLIST command shows your MAPLIST.map file; and the OTHER MAPS command"
  57.    PRINT " shows all MAP borders so you can see how your new map fits in.  You may use F3"
  58.    PRINT " and F4 keys to select smaller or larger map borders to draw."
  59.    PRINT
  60.    LOCATE 25, 1: PRINT " HIT ANY KEY to continue...";
  61.    GOSUB GetChar: BEEP
  62.   
  63. Info3:   COLOR 15, 3: CLS
  64.    PRINT " PAGE 3 INSTRUCTIONS:  Using GPS Track History Files to draw maps!"
  65.    PRINT
  66.    PRINT " To aid in creating accurate maps, W7KKE in California added routines to MAPFIX"
  67.    PRINT " so that you can overlay a Track History file onto the map you are constructing."
  68.    PRINT " This is an excellent tool for correcting your maps to real GPS data."
  69.    PRINT ""
  70.    PRINT " After you have loaded your map, type alt-G (GPS Track History') and enter the"
  71.    PRINT " history filename.  This will overlay the track history file.  You may then"
  72.    PRINT " use the normal MAPFIX.bas routines to move map segments and add so that the"
  73.    PRINT " map will agree with the GPS data contained in the track history file.  This"
  74.    PRINT " is especially useful with the history files saved by a laptop during mobile"
  75.    PRINT " GPS operations."
  76.    PRINT
  77.    PRINT " CAUTION:  Since GPS data is only accurate to 100 yards due to the effects of"
  78.    PRINT " Selective Availability, I would avoid using GPS data explicitely below about"
  79.    PRINT " the 2 mile range.  For this reason, I make the size of the GPS positions "
  80.    PRINT " expand below the 2 mile range to roughly approximate the size of the 100 yard"
  81.    PRINT " error circle."
  82.    PRINT
  83.    PRINT " Also note that you can now START a NEW map from scratch, without using the"
  84.    PRINT " text editor to enter the initial configuration data.  Just type NEW instread "
  85.    PRINT " of a MAPfilename when starting up the program."
  86.    LOCATE 25, 1: PRINT " HIT ANY KEY to continue...";
  87.    GOSUB GetChar: BEEP
  88.  
  89. Info4:   COLOR 15, 7: CLS
  90.    PRINT " PAGE 4 INSTRUCTIONS FOR USING A DIGITIZER:": PRINT
  91.    PRINT " In version 3.07B, MAPFIX.bas can now interface to a digitizer tablet or table"
  92.    PRINT " so that maps can be drawn directly from paper maps.  These routines were"
  93.    PRINT " developed by W7KKE and added to MAPFIX.bas in December 1993."
  94.    PRINT
  95.    PRINT " To use a digitizer, first you must hit the ALT-O command to open the COMM"
  96.    PRINT " PORT for the digitizer.  This command also lets you test the digitizer while"
  97.    PRINT " testing the alignment of the map on the digitizer surface.  It then prompts"
  98.    PRINT " you to identify the upper left and lower right corners of the map, in order"
  99.    PRINT " to calibrate the digitizer to the latitude, longitude and scale of the map."
  100.    PRINT
  101.    PRINT " From this point on, the button on the digitizer mouse is almost identical to"
  102.    PRINT " the ALT-A command for ADDing a point.  To start a new map feature, however,"
  103.    PRINT " for the digitizer, you DO NOT use the ALT-N NEW command, but you should use"
  104.    PRINT " the ALT-B BEGIN command.  For More information, see the README.DIG file."
  105.    PRINT
  106.    PRINT " To speed up the map drawing during editing, I no longer erase and re-draw"
  107.    PRINT " the entire map with each new point.  I simply draw just the new line segment."
  108.    PRINT " Sometimes, especially when you move, or add a line, this leaves an old line"
  109.    PRINT " segment, where there actually is no longer a line.  You can always celan up"
  110.    PRINT " the map by just hitting the space bar to force a new map..."
  111.    PRINT
  112.  
  113.   
  114.    Display$ = "UNKnown"
  115.    RETURN
  116.  
  117. GetChar: a$ = "": DO UNTIL a$ <> "": a$ = INKEY$: LOOP: RETURN
  118.  
  119. BEGIN: GOSUB Info:
  120.    PRINT " HIT ANY KEY to proceed onto the HELP screen...";
  121.    GOSUB GetChar
  122.  
  123.    DIM x%(5 * MaxNumPoints), y%(5 * MaxNumPoints)
  124.    REM MAP coordinates               **** THESE ARE BIGGER THAN APRS ***
  125.    DIM LN$(MaxNumLines) ' (no limit in APRS)  **** SO YOU CAN MANIPULATE BIG MAPS
  126.    nn = 2 * MaxNumLABELS
  127.    DIM ML$(nn), MLa(nn), MLo(nn), MLr(nn) 'Map Labels, lengths and coordinates
  128.    nn = 2 * MaxNumMAPS
  129.    DIM MapName$(nn), LATcen(nn), LONcen(nn), MapMax(nn), Comment$(nn)
  130.    RdsOn = -1: Labls = -1: Tags = -1: KP = 1: Changed = 0: MapSize = 256
  131.    i = 1000
  132.    DIM HLAT(i), HLONG(i)'For lat/longs from big GPS history files
  133.  
  134. INIT: ON ERROR GOTO ErrorTrap
  135.    ScrnType$ = "EGA": Ycen = 200: Yfactr = 1: YfacTXT = 350 / 350: SCREEN 9
  136.    IF ScrnType$ = "EGA" THEN COLOR 15, 0
  137.    REM ScrnType$ = "CGA": Yfactr=200/400:Ycen = 200*Yfactr: SCREEN 2
  138.   
  139. Display$ = "HELP": GOSUB HELP: GOSUB LoadMap
  140. REM ON ERROR GOTO 0
  141.  
  142. Main: GOSUB DrwMPaCur
  143.   DO
  144. GoAgain: Fault = 0
  145.     IF Digitizer THEN
  146.          IF LOC(1) > 9 THEN
  147.             GOSUB GetXY: GOSUB Cursor
  148.             IF Btn <> 3 THEN GOSUB AddPoint
  149.          END IF
  150.     END IF
  151.     a$ = INKEY$
  152.     IF a$ <> "" THEN
  153.         a$ = UCASE$(a$): Key$ = a$
  154.         IF a$ = "S" THEN GOSUB labels
  155.         IF a$ = "L" THEN Labls = NOT Labls
  156.         IF a$ = "T" THEN Tags = NOT Tags
  157.         IF a$ = "F" THEN GOSUB FindPoint
  158.         IF a$ = "G" THEN GOSUB CurToPoint: GOSUB CurDrwMap
  159.         IF a$ = "H" THEN
  160.            IF Display$ <> "HELP" THEN
  161.               GOSUB HELP
  162.            ELSE GOSUB Info
  163.               LOCATE 25, 1: PRINT " H for HELP or SPACE BAR for map..."; : a$ = ""
  164.            END IF
  165.         END IF
  166.         IF a$ = "B" THEN GOSUB BoxPPD
  167.         IF a$ = "U" THEN GOSUB GetUSGS
  168.         IF a$ = "D" THEN GOSUB MapDIR
  169.         IF a$ = "M" THEN GOSUB ListMAPlist
  170.         IF a$ = "O" THEN GOSUB DrwAndShow
  171.         IF a$ = "N" THEN GOSUB NextLine: GOSUB Cursor
  172.         IF a$ = "P" THEN GOSUB Previous: GOSUB Cursor
  173.         IF a$ = "Q" THEN GOSUB QUIT
  174.         IF a$ = "R" THEN Z = 2: LNptr = 1
  175.         IF a$ = "T" THEN GOSUB Scrunch
  176.         IF a$ = " " THEN Display$ = "MAP": GOSUB DrwMPaCur
  177.         IF a$ = "+" THEN Z = Z + 1: GOSUB MapPoint ' moves to next map point
  178.         IF a$ = "-" THEN Z = Z - 1: GOSUB MapPoint ' moves backwards
  179.         IF a$ = CHR$(18) THEN USGS = NOT USGS: GOSUB ReDraw
  180.        
  181.         B$ = "": IF LEN(a$) = 2 THEN B$ = RIGHT$(a$, 1): REM process arrow & special keys
  182.         IF B$ = "I" THEN RS = RS * 2: GOSUB CurDrwMap: REM change scale
  183.         IF B$ = "Q" THEN RS = RS / 2: GOSUB CurDrwMap
  184.         IF B$ = CHR$(132) THEN RS = RS * 8: GOSUB CurDrwMap: REM change scale by factor of 4
  185.         IF B$ = "V" THEN RS = RS / 8: GOSUB CurDrwMap
  186.         IF B$ = "G" THEN GOSUB CurDrwMap 'Home key
  187.         IF a$ = "7" THEN CDX = LONo: CDY = LATo: GOSUB DrwMPaCur 'ShiftHOME
  188.         IF B$ = "O" THEN CDX = LONcen: CDY = LATcen: GOSUB DrwMPaCur 'End Key
  189.         IF B$ = "M" THEN CPX = CPX - 4 / (Sfac): GOSUB Cursor
  190.         IF B$ = "K" THEN CPX = CPX + 4 / (Sfac): GOSUB Cursor
  191.         IF B$ = "H" THEN CPY = CPY + 4 / (Sfac): GOSUB Cursor
  192.         IF B$ = "P" THEN CPY = CPY - 4 / (Sfac): GOSUB Cursor
  193.         REM Here are the special MapFIx routines
  194.         IF B$ = CHR$(30) THEN GOSUB AddPoint               'alt-ADD point
  195.         IF B$ = CHR$(48) AND Digitizer THEN GOSUB NewFeature'alt-BEGIN
  196.         IF B$ = CHR$(34) THEN GOSUB LoadHst                'alt-GPS hstry file
  197.         IF B$ = CHR$(50) THEN GOSUB MakePT: IF NOT USGS THEN GOSUB DrawMap 'MOVE point to cursor
  198.         IF B$ = CHR$(32) THEN GOSUB DelPT                  'alt-DELete point
  199.         IF B$ = CHR$(38) THEN GOSUB AddLabel               'alt-ADD LABEL
  200.         IF B$ = CHR$(46) THEN GOSUB NewCenter              'alt-CENTER
  201.         IF B$ = CHR$(36) THEN GOSUB Join                   'alt-JOIN
  202.         IF B$ = CHR$(37) THEN GOSUB KillF                  'alt-KILL Feature
  203.         IF B$ = CHR$(19) THEN GOSUB MapRange               'alt-RANGE
  204.         IF B$ = CHR$(20) THEN GOSUB TRIM                   'alt-TRIM
  205.         IF B$ = CHR$(49) THEN GOSUB NewFeature             'alt-NEW Feature
  206.         IF B$ = CHR$(24) THEN GOSUB DigiInit: GOSUB DrawMap'alt-OPEN dgtzr COM
  207.         IF B$ = CHR$(31) THEN GOSUB Scrunch                'alt-SCRUNCH
  208.         IF B$ = CHR$(22) THEN GOSUB GetUSGS                'alt-U
  209.         IF B$ = CHR$(61) THEN                              'F3 for smaller Maps
  210.            MapSize = MapSize / 2: IF MapSize < 1 THEN MapSize = 1
  211.            GOSUB ShowMaps
  212.            END IF
  213.         IF B$ = CHR$(62) THEN                              'F4 for larger Maps
  214.            MapSize = MapSize * 2: IF MapSize > 1000 THEN MapSize = 1000
  215.            GOSUB DrwAndShow
  216.         END IF
  217.         IF a$ = CHR$(19) THEN GOSUB SaveMap
  218.            
  219.         IF a$ = "6" THEN CPX = CPX - 20 / (Sfac): GOSUB Cursor'SHIFT Cursor by 4
  220.         IF a$ = "4" THEN CPX = CPX + 20 / (Sfac): GOSUB Cursor
  221.         IF a$ = "8" THEN CPY = CPY + 20 / (Sfac): GOSUB Cursor
  222.         IF a$ = "2" THEN CPY = CPY - 20 / (Sfac): GOSUB Cursor
  223.         
  224.     END IF
  225.   LOOP
  226.   SYSTEM 'you should never get here
  227.  
  228. ReDraw: LOCATE 1, 30
  229.         IF USGS THEN PRINT "NO REDRAW     ":  ELSE PRINT "REDRAW ENABLED"
  230.         RETURN
  231.  
  232. QUIT: a$ = "Y"
  233.       IF Changed THEN
  234.          GOSUB BoxLine23
  235.          PRINT "**** MAP HAS BEEN MODIFIED"; Changed; "TIMES BUT NOT SAVED!!!  SAVE NOW? (Y)";
  236.          INPUT a$
  237.       IF UCASE$(a$) <> "N" THEN GOSUB SaveMap
  238.       END IF
  239.       SYSTEM
  240.  
  241. TRIM: GOSUB BoxLine23
  242.       CLS : PRINT "TRIM ALL POINTS AND LABELS OUTSIDE OF MAPRANGE"
  243.       PRINT
  244.       PRINT "This command will remove all points and labels that are outside of the white"
  245.       PRINT "map border.  You can change the location of this map border by using"
  246.       PRINT "the CENTER command (alt-C) and by changing the RANGE using alt-R."
  247.       PRINT : PRINT
  248.       PRINT "No map feature will be completely eliminated..."
  249.       PRINT
  250.       PRINT "The first and last point of any FEATURE will be retained, so the"
  251.       PRINT "result will be long single lines for all FEATURES outside the map border."
  252.       PRINT "Use the KILL FEATURE (alt-K) to eliminate those lines and use the MOVE"
  253.       PRINT "command (alt-M) to move any far away points closer to the border."
  254.       PRINT : PRINT
  255.       PRINT "You might consider stopping now and doing a SAVE (ctrl-S) before proceeding."
  256.       PRINT
  257.       PRINT "ALSO, THIS DOES NOT WORK FOR POINTS WITH NEGATIVE VALUES!  Be sure  your"
  258.       PRINT "selected area is below and to right of ORIGIN.  If not, run MAPCNVRT.bas."
  259.       PRINT : PRINT
  260.       INPUT "Are you ready to proceed? (Y/N) (N)"; ans$
  261.       GOSUB DrawMap
  262.       IF UCASE$(ans$) <> "Y" THEN RETURN
  263.       C = 0: LOCATE 23, 1: PRINT "Processing...";
  264.       REM dx and dy are num pix of center of map
  265.       REM bx and by are borders of map based on MapRng
  266.       by = ppdV * MapRng / 60
  267.       bx = by / Lfac
  268.       FOR Z = 1 TO nmp - 4
  269.          IF x%(Z) = 0 THEN Z = Z + 2
  270.          IF x%(Z) > dx + bx OR y%(Z) > dy + by THEN bad = 1 ELSE bad = 0
  271.          IF x%(Z) < dx - bx OR y%(Z) < dy - by THEN bad = 1
  272.          IF bad AND x%(Z - 1) <> 0 AND x%(Z + 1) <> 0 THEN
  273.             GOSUB DelPT: Z = Z - 1
  274.             C = C + 1
  275.          END IF
  276.       NEXT Z
  277.       LOCATE 23, 1: PRINT "Now removing labels...";
  278.       FOR i = 1 TO nml: REM now eliminate all labels outside
  279.           bad = 0: Xm = MapRng / (60 * Lfac): Ym = MapRng / 60
  280.           IF MLo(i) > LONcen + Xm OR MLa(i) > LATcen + Ym THEN bad = 1
  281.           IF MLo(i) < LONcen - Xm OR MLa(i) < LATcen - Ym THEN bad = 1
  282.           IF bad = 1 THEN
  283.              FOR j = i TO nml
  284.                  ML$(j) = ML$(j + 1): MLa(j) = MLa(j + 1)
  285.                  MLo(j) = MLo(j + 1): MLr(j) = MLr(j + 1)
  286.              NEXT j: nml = nml - 1: PRINT ".";
  287.           END IF
  288.       NEXT i
  289.       GOTO DrawMap
  290.  
  291.  
  292. FindPoint: CurX = INT(.5 + dx + (CUX - 320) / KP)
  293.            CurY = INT(.5 + dy + (CUY - Ycen) / KP)
  294.     GOSUB BoxLine23: PRINT "SEARCHING THROUGH ALL POINTS IN FILE...";
  295.     FOR j = 0 TO 30             ' Go through abt 20 times lookin pt.
  296.         IF j > 10 THEN j = j + 1' first with 0 delta, then bigger
  297.         LNctr = 0: PRINT ".";
  298.         FOR i = 1 TO nmp
  299.             IF x%(i) = 0 THEN LNctr = LNctr + 1
  300.             IF x%(i) > CurX - j AND x%(i) < CurX + j THEN
  301.                IF y%(i) > CurY - j AND y%(i) < CurY + j THEN
  302.                   Z = i: LNptr = LNctr: GOSUB CurToPoint
  303.                   j = 99: i = nmp
  304.                END IF
  305.             END IF
  306.         NEXT i:
  307.     NEXT j
  308.     IF j < 99 THEN PRINT "None found!": RETURN
  309.     GOSUB MapPoint: SavClr = 0: RETURN
  310.                   
  311. NewFeature: LOCATE 24, 1: PRINT SPACE$(27); : GOSUB BoxLine23
  312.             INPUT "Enter reference name for new feature"; a$
  313.             IF a$ = "" THEN RETURN
  314.             LOCATE 25, 1
  315.             FOR i = 0 TO 14
  316.             PRINT RIGHT$(" " + MID$(STR$(i + 1), 2), 2); "   ";
  317.             LINE (16 + i * 40, 335 * YfacTXT)-(40 + i * 40, 349 * YfacTXT), i + 1, BF
  318.             NEXT i
  319.             GOSUB BoxLine23
  320.             INPUT "Select color (4,7,10-Hwys 11-Water 12-Hwy 13-Spcl 14-City)"; B$
  321.             SavClr = VAL(B$): IF SavClr > 15 OR SavClr < 1 THEN RETURN
  322.             GOSUB BeginF
  323.             GOSUB BoxLine23: LOCATE 25, 1: PRINT SPACE$(80); : LOCATE 25, 1
  324.             IF RIGHT$(Key$, 1) = CHR$(48) THEN
  325.                PRINT "NOW USE DIGITIZER TO ADD NEW POINTS TO THIS FEATURE...";
  326.                GOSUB GetXY: GOSUB Cursor
  327.             ELSE
  328.                PRINT "NOW MOVE CURSOR AND USE ALT-A TO ADD POINTS TO THIS NEW FEATURE...";
  329.             END IF
  330.             GOSUB MakePT
  331.             RETURN
  332.            
  333. BeginF: x%(nmp) = 0: y%(nmp) = SavClr   'Store feature color 0,c
  334.         LN$(LNi + 1) = LN$(LNi): LNptr = LNi'Bump up present LN$ comment
  335.         LN$(LNi) = a$: LNi = LNi + 1'Store feature name
  336.         nmp = nmp + 1: Z = nmp
  337.         nmp = nmp + 1: x%(nmp) = 0: y%(nmp) = 0'nmp points to ending 0,0
  338.         RETURN
  339.  
  340. CanclF: nmp = nmp - 2: Z = Kz
  341.         LNi = LNi - 1: LN$(LNi) = LN$(LNi + 1): RETURN
  342.  
  343. NewCenter: LATcen = CPY: LONcen = CPX: Changed = Changed + 1: GOTO CurDrwMap
  344.  
  345. MapRange: GOSUB BoxLine23: INPUT "Enter map range"; a$
  346.           IF VAL(a$) <> 0 THEN MapRng = VAL(a$)
  347.           Changed = Changed + 1: GOTO DrwMPaCur
  348.  
  349. AddPoint: nmp = nmp + 1: Z = Z + 1
  350.           FOR i = nmp TO Z STEP -1
  351.               x%(i) = x%(i - 1): y%(i) = y%(i - 1)
  352.           NEXT
  353.           GOSUB MakePT
  354.           IF SavClr = 0 THEN GOTO DrawMap
  355.           s = Z - 1: LineColor = SavClr: GOTO DP
  356.  
  357. MakePT: x%(Z) = dx + (CUX - 320) / (KP * Hfac)
  358.         y%(Z) = dy + (CUY - Ycen) / KP
  359.         Changed = Changed + 1
  360.         GOTO MapPoint
  361.  
  362. CurToPoint:
  363.      CPX = CDX - (x%(Z) - dx) / ppdV
  364.      CPY = CDY - (y%(Z) - dy) / (ppdV * Yfactr)
  365.      GOTO Cursor
  366.  
  367. DelPT: GOSUB DelZ
  368.        REM if 1st pt, it stays as 1st pt
  369.  
  370.        IF x%(Z) = 0 THEN Z = Z - 1: REM if end pt, it stays as end
  371.        IF x%(Z + 1) = 0 AND x%(Z - 1) = 0 THEN 'It is LAST point
  372.           GOSUB Kline: LNptr = LNptr - 1       'So Kill Line
  373.           GOSUB DelZ                           'And Kiil it
  374.           Z = Z - 1: GOSUB DelZ: Z = Z - 1     'Kill 0,color
  375.        END IF                                  'and -1 to end point
  376.        IF B$ = CHR$(32) AND USGS = 0 THEN GOSUB DrawMap ELSE GOSUB MapPoint
  377.        RETURN
  378.  
  379. DelZ: nmp = nmp - 1
  380.       FOR i = Z TO nmp
  381.           x%(i) = x%(i + 1): y%(i) = y%(i + 1)
  382.       NEXT: Changed = Changed + 1: RETURN
  383.  
  384. NextLine: IF Z >= nmp - 1 THEN Z = nmp - 1: BEEP: RETURN
  385.           DO UNTIL x%(Z) = 0: Z = Z + 1: LOOP
  386.           IF Z < nmp - 1 THEN Z = Z + 1: LNptr = LNptr + 1
  387.           SavClr = 0: GOTO MapPoint
  388. Previous: DO UNTIL Z = 1 OR x%(Z) = 0: Z = Z - 1: LOOP
  389.           IF Z > 3 THEN Z = Z - 1: LNptr = LNptr - 1
  390.           SavClr = 0: GOTO MapPoint
  391.  
  392. KillF: Bi = Z: Changed = Changed + 1
  393.        DO UNTIL x%(Bi) = 0: Bi = Bi - 1: LOOP: Z = Bi + 1
  394.        REM Stop at Beginning (0) point of the feature to kill
  395.        ni = Bi + 1' Now scan for next feature
  396.        DO UNTIL x%(ni) = 0: ni = ni + 1: LOOP
  397.        REM now move down rest of array to fill
  398.        DO UNTIL ni = nmp + 1
  399.            x%(Bi) = x%(ni): y%(Bi) = y%(ni)
  400.            Bi = Bi + 1: ni = ni + 1
  401.        LOOP
  402.        nmp = nmp - (ni - Bi): y%(nmp) = 0
  403.        GOSUB Kline
  404.        GOTO DrawMap
  405.      
  406. Kline: FOR i = LNptr TO LNi
  407.            LN$(i) = LN$(i + 1)
  408.            NEXT i
  409.            LNi = LNi - 1
  410.            RETURN
  411.  
  412. MapPoint:
  413.      IF Z < 2 THEN Z = 2: LNptr = 1: BEEP: SavClr = 0
  414.      IF Z > nmp - 1 THEN Z = Z - 1: BEEP: SavClr = 0
  415.      IF x%(Z) = 0 THEN
  416.         IF a$ = "-" THEN
  417.              LNptr = LNptr - 1: Z = Z - 1
  418.         ELSE LNptr = LNptr + 1: Z = Z + 1
  419.         END IF: SavClr = 0
  420.      END IF
  421.      IF LNptr < 0 THEN LNptr = 0
  422.      IF Display$ = "MAP" THEN
  423.           LOCATE 22, 1
  424.           PRINT "Fture#"; LNptr; TAB(12); LEFT$(LN$(LNptr) + "            ", 12);
  425.      END IF
  426. DrwMpPt: IF Display$ <> "MAP" THEN RETURN
  427.      CIRCLE (Xtest, Ytest), 10, 0 'Erase old circle
  428.      Xtest = 320 + KP * (x%(Z) - dx) * Hfac
  429.      Ytest = Ycen + KP * (y%(Z) - dy) * Yfactr
  430.      CIRCLE (Xtest, Ytest), 10, 15
  431.      
  432.      LOCATE 23, 1: PRINT "MapPt#"; Z;
  433.      IF Z > 999 THEN PRINT TAB(13); "val:";  ELSE PRINT TAB(12); "vals:";
  434.      PRINT TAB(17); x%(Z); TAB(23); y%(Z)
  435.      RETURN
  436.  
  437. AddLabel: nml = nml + 1
  438.           MLa(nml) = CPY: MLo(nml) = CPX
  439.           GOSUB BoxLine23: INPUT "Enter Label Name"; a$: ML$(nml) = a$
  440.           GOSUB BoxLine23: INPUT "Begin displaying label at what range?"; a$
  441.           a = VAL(a$): IF a <> 0 THEN MLr(nml) = a:  ELSE MLr(nml) = 2048
  442.           Changed = Changed + 1: GOTO labels
  443.  
  444. BoxLine23: LOCATE 23, 1: PRINT SPACE$(80); : LOCATE 23, 1: RETURN
  445.  
  446. ErrorTrap: Fault = ERR: 'Error handling routine
  447.            IF ERR = 57 THEN PRINT "  I/O-error-User-logoff"; : RESUME
  448.            IF ERR = 69 THEN PRINT "  Comm-buffer-overflow"; : RESUME
  449.            IF ERR = 53 THEN PRINT "  file-"; F$; "-not-found": CLOSE : RESUME NEXT
  450.            IF ERR = 62 THEN RESUME NEXT
  451.            IF ERR = 52 THEN RESUME NEXT
  452.            IF ERR = 55 THEN RESUME NEXT
  453.            IF ERR = 2 THEN PRINT "SYNTAX-error"
  454.            IF ERR = 70 THEN PRINT " WRITE PROTECTED!...": RESUME NEXT
  455.            IF ERR = 76 THEN PRINT "Wrong Path!": RESUME NEXT
  456.            RESET
  457.            PRINT : PRINT "Error beyond repair. Number = "; ERR;
  458.            INPUT "Hit RETURN to return to DOS"; a$
  459.            SYSTEM
  460.  
  461. MapDIR: CLS : PRINT "MAP FILES DIRECTORY": PRINT
  462.          PRINT "To display MAP files, please enter the path to your xxxxxxx.MAP files."
  463.          PRINT "For example, the default '\APRS\MAPS\*.MAP' will show all maps in the APRS"
  464.          PRINT "directory.  Similarly '*.map' will search your present QB directory."
  465.          PRINT "For any other path, enter the full file specification.": PRINT
  466.          F$ = "\aprs\MAPS\*.map"
  467.          PRINT "Enter Filespec for searching the DIRECTORY ("; F$; ")";
  468.          INPUT a$: IF a$ <> "" THEN F$ = a$
  469.          PRINT : PRINT : FILES F$
  470.          RETURN
  471.  
  472.  
  473. LoadMap: 'Maps are drawn to the default EGA resolution of 640 x 400 (350)
  474. Again: GOSUB BoxLine23
  475.        INPUT " Enter map FILENAME, or NEW, or ? for a list, or Q to quit)"; a$
  476.        a$ = UCASE$(a$): IF a$ = "" THEN GOTO Again
  477.        IF a$ = "Q" THEN SYSTEM
  478.        IF a$ = "?" THEN GOSUB MapDIR: GOTO Again
  479.        IF a$ = "NEW" THEN Key$ = "NEW": GOSUB NewMap: RETURN
  480.        a = INSTR(3, a$, "."): IF a = 0 THEN a$ = a$ + ".MAP"
  481.        MapFile$ = a$: F$ = MapFile$: OPEN F$ FOR INPUT AS #3
  482.        IF Fault = 53 THEN Fault = 0: PRINT : CLOSE #3: GOTO Again
  483.        GOSUB BoxLine23: PRINT " Loading "; F$; "..."
  484.        INPUT #3, LATo: LINE INPUT #3, LATtext$
  485.        INPUT #3, LONo: LINE INPUT #3, LONtext$
  486.        INPUT #3, ppdV: LINE INPUT #3, VS$'Pixels per degree horiz
  487.        INPUT #3, LATcen: LINE INPUT #3, LATcen$
  488.        INPUT #3, LONcen: LINE INPUT #3, LONcen$
  489.        INPUT #3, MapRng: LINE INPUT #3, MapRng$
  490.        INPUT #3, MinRng: LINE INPUT #3, MR$
  491.        LINE INPUT #3, TextLine$ ' Line of comments or instrutcitons
  492.        IF RIGHT$(TextLine$, 6) = "CD ROM" THEN USGS = -1
  493.        RS = 2 ^ INT(LOG(MapRng) / LOG(2))'Rng is intgr of VERTrng
  494.        i = 0: LNi = 0:
  495.      
  496.      DO WHILE NOT EOF(3)
  497.         i = i + 1: INPUT #3, x%(i), y: y%(i) = y * Yfactr
  498.         IF x%(i) = 0 AND NOT EOF(3) THEN ' Get line color & store with x=0
  499.            INPUT #3, y%(i): LNi = LNi + 1: LINE INPUT #3, LN$(LNi)' Save line name
  500.            IF y = -1 THEN GOSUB LoadLabels ' All labels listed at end of file
  501.            END IF
  502.         LOOP: nmp = i  'nmp points to 0,-1 that ends all data (but the value
  503.                        'of X% and y% are 0,0 until file is saved.
  504.      LET CDY = LATcen: CDX = LONcen'Center display on ORIGIN
  505.      LET CPX = CDX: CPY = CDY 'Cursor Posn to Center of Display
  506.      LET Z = 2: LNptr = 1: REM start at first point and first line segment
  507.      CLOSE #3: RETURN:     REM first X% value is map color.  2nd val is 1st pt
  508.  
  509.  
  510.  
  511. LoadLabels: k = 0
  512.      DO WHILE NOT EOF(3)
  513.         k = k + 1: INPUT #3, ML$(k), MLa(k), MLo(k), MLr(k)
  514.      LOOP
  515.      IF MLa(k) = 0 OR MLo(k) = 0 THEN nml = k - 1 ELSE nml = k
  516.      RETURN
  517.        
  518. SaveMap: GOSUB BoxLine23
  519.    PRINT "Enter file name to save if other than "; MapFile$;
  520.    INPUT a$: IF a$ <> "" THEN MapFile$ = a$
  521.    F$ = MapFile$
  522.    GOSUB BoxLine23: PRINT "Saving map to file named "; F$; " ..."
  523.    OPEN F$ FOR OUTPUT AS #4
  524.    IF Fault = 70 THEN CLOSE #4: GOTO SaveMap
  525.    PRINT #4, LATo; ","; LATtext$
  526.    PRINT #4, LONo; ","; LONtext$
  527.    PRINT #4, ppdV; ","; VS$
  528.    PRINT #4, LATcen; ","; LATcen$
  529.    PRINT #4, LONcen; ","; LONcen$
  530.    PRINT #4, MapRng; ","; MapRng$
  531.    PRINT #4, MinRng; ","; MR$
  532.    PRINT #4, TextLine$
  533.    j = 1
  534.    FOR i = 1 TO nmp
  535.        IF x%(i) <> 0 THEN WRITE #4, x%(i), INT((y%(i) / Yfactr) + .5)
  536.        IF x%(i) = 0 AND i = nmp THEN PRINT #4, " 0,-1"
  537.        IF x%(i) = 0 AND i <> nmp THEN
  538.           PRINT #4, "0,0"
  539.           PRINT #4, y%(i); ","; LN$(j): j = j + 1
  540.        END IF
  541.    NEXT i
  542.    PRINT #4, "0,"; LN$(LNi)
  543.    FOR k = 1 TO nml
  544.        PRINT #4, ML$(k); ","; : WRITE #4, MLa(k), MLo(k), MLr(k)
  545.    NEXT k: CLOSE #4: LOCATE 24, 1:
  546.    Changed = 0
  547.    IF nmp > MaxNumPoints OR nml > MaxNumLABELS THEN
  548.       CLS : LOCATE 9, 29: PRINT "CAUTION!": PRINT : PRINT
  549.       IF nmp > MaxNumPoints THEN
  550.          PRINT "            The number of points,"; nmp; "is greater than"; MaxNumPoints
  551.       END IF
  552.       IF nml > MaxNumLABELS THEN
  553.          PRINT "            The number of LABELS,"; nml; "is greater than"; MaxNumLABELS
  554.       END IF
  555.       LOCATE 18, 12
  556.       PRINT " Therefore this map will not work with APRS (yet) "
  557.       LOCATE 23, 1: INPUT "HIT Enter to continue..."; a$
  558.    END IF: GOTO DrwMPaCur
  559.  
  560. CurDrwMap: CDX = CPX: CDY = CPY: GOTO DrawMap: REM Re-center at CURSOR location
  561.  
  562. DrwMPaCur: CPX = CDX: CPY = CDY: GOSUB DrawMap
  563.            REM After drawing map, Put cursor at center
  564.            RETURN
  565.  
  566. DrawMap: Display$ = "MAP": COLOR 15, 0
  567.    'Draw to range scale RS and center display CDX and CDY
  568.    'Original Map was 40 pix-per-deg Horiz and 20 vert for 200 display
  569.    'Now ppdH and ppdV are variables.  The scaling factor KP is 1 for
  570.    'the original map scale.
  571.    DO WHILE RS < 320 / ppdV: RS = RS * 2: LOOP
  572.    IF RS > 8192 THEN RS = 8192
  573.    KP = 100 * 100 / (RS * ppdV)'This is to scale it down from the 120 maps
  574.    Sfac = 50 * 200 / RS        'Till 307 had been 100*120
  575.  
  576.    Lfac = COS(CDY / 57.296)
  577.    Hfac = (640 / 350) * (3 / 4) * Lfac
  578.    dx = ppdV * (LONo - CDX)
  579.    dy = ppdV * (LATo - CDY)
  580.         
  581.    CLS : LOCATE 1, 2: PRINT "Redrawing Map"
  582.    REM first put ORIGIN and map CENTER on the map
  583.    LINE (320 - KP * dx, Ycen - KP * dy)-(960 - KP * dx, Ycen - KP * dy), 14
  584.    LINE (320 - KP * dx, Ycen - KP * dy)-(320 - KP * dx, 3 * Ycen - KP * dy), 14
  585.    CMX = 320 + Sfac * (CDX - LONcen) * Hfac'new
  586.    CMY = Ycen + Sfac * (CDY - LATcen) * Yfactr
  587.    LINE (CMX - 27, CMY)-(CMX + 27, CMY), 14
  588.    LINE (CMX, CMY - 20)-(CMX, CMY + 20), 14
  589.    CIRCLE (CMX, CMY), 10, 14
  590.    CIRCLE (320 - KP * dx, Ycen - KP * dy), 12, 14
  591.    s = 0: GOSUB MapPoint: REM Redraw MapPoint
  592.    StrtPt = -1
  593. DP: FOR i = s TO nmp - 1
  594.        x = 320 + KP * (x%(i) - dx) * Hfac
  595.        y = Ycen + KP * (y%(i) - dy) * Yfactr
  596.        X1 = 320 + KP * (x%(i + 1) - dx) * Hfac
  597.        Y1 = Ycen + KP * (y%(i + 1) - dy) * Yfactr
  598.        IF StrtPt = -1 THEN CIRCLE (x, y), 3, 9: StrtPt = 0
  599.        IF x%(i + 1) <> 0 THEN
  600.           IF RdsOn OR LineColor <> 12 THEN LINE (x, y)-(X1, Y1), LineColor
  601.           IF i = Z THEN SavClr = LineColor
  602.        ELSE
  603.           LINE (x - 3, y - 3)-(x + 3, y + 3), 10, B: StrtPt = -1
  604.           LineColor = y%(i + 1): i = i + 1
  605.           IF Display$ = "SHOW" AND LineColor > 8 THEN LineColor = LineColor - 8
  606.        END IF
  607.    NEXT i
  608.    GOSUB Cursor
  609.    GOSUB ReDraw
  610.    REM MapPoint went here
  611.    GOSUB DrawHist: REM draw GPS history track
  612.    IF Display$ = "SHOW" THEN
  613.       GOSUB ShowMaps
  614.    ELSE
  615.       LOCATE 25, 1: PRINT "Use +/- to move MAPpoint.  N/P for Next/Previous Feature.  H for HELP!.";
  616.       LOCATE 1, 61
  617.       PRINT "POINTS"; nmp; "= "; INT((nmp / MaxNumPoints) * 100); "%";
  618.       LOCATE 2, 61
  619.       PRINT "LABELS "; nml; "= "; INT((nml / MaxNumLABELS) * 100); "%";
  620.       LOCATE 3, 61: PRINT "CENTER  "; MID$(STR$(LATcen), 2, 5);
  621.       LOCATE 3, 75: PRINT MID$(STR$(LONcen), 2, 5)
  622.       LOCATE 4, 61: PRINT "SCALE   (ppd)"; ppdV
  623.       LOCATE 5, 69: PRINT "Range"; LEFT$(STR$(MapRng), 5)
  624.    END IF
  625.  
  626. labels:
  627.    IF Labls THEN
  628.       FOR i = 1 TO nml ' Now plot labels on map
  629.       IF RS <= MLr(i) OR Key$ = "S" THEN
  630.          LET x = 320 + Sfac * (CDX - MLo(i)) * Hfac'new
  631.          LET y = Ycen + Sfac * (CDY - MLa(i)) * Yfactr
  632.          IF Tags AND y > 14 * Yfactr AND y < 350 * Yfactr AND x > 8 * (LEN(ML$(i)) + 1) AND x < 632 THEN
  633.             LOCATE y / (14 * Yfactr), (x / 8) - LEN(ML$(i)): PRINT ML$(i);
  634.             END IF
  635.          END IF
  636.       NEXT i
  637.    END IF
  638.    GOSUB ShowMap: RETURN
  639.          
  640. ShowMap: REM this shows the map boarder of the loaded map
  641.     x = 320 + KP * (CDX - LONcen) * ppdV * Hfac'new
  642.     y = Ycen + KP * (CDY - LATcen) * ppdV * Yfactr
  643.       by = MapRng * Sfac * Yfactr / 60
  644.       bx = by * 640 / (400 * Yfactr) * Lfac'old
  645.       C = 15
  646.     LINE (x - bx, y - by)-(x + bx, y + by), C, B
  647.     RETURN
  648.  
  649. Cursor: CIRCLE (CUX, CUY), 4, 0
  650.      CUX = 320 + Sfac * (CDX - CPX) * Hfac'new
  651.      CUY = Ycen + Sfac * (CDY - CPY) * Yfactr
  652.      CIRCLE (CUX, CUY), 4, 14
  653.      x = INT(CPX): y = INT(CPY): Xm = (CPX - x) * 60: Ym = (CPY - y) * 60
  654.      x$ = RIGHT$(STR$(x), 3) + " "
  655.      LOCATE 1, 2: PRINT "RNG"; RIGHT$("   " + STR$(RS), 4) + "      "
  656.      LOCATE 2, 2: PRINT "LAT "; y; MID$(STR$(Ym) + "   ", 2, 5)
  657.      LOCATE 3, 2: PRINT "LON "; x$; MID$(STR$(Xm) + "   ", 2, 5)
  658.     
  659.      LOCATE 24, 1: PRINT "Cursor coordnts:"; TAB(17);
  660.      PRINT INT(.5 + dx + (CUX - 320) / KP); TAB(23); INT(.5 + dy + (CUY - Ycen) / KP);
  661.      REM LOCATE 24, 55: PRINT "Degrees: ";
  662.      REM PRINT LEFT$(STR$(CPY) + " ", 7); LEFT$(STR$(CPX) + "   ", 7);
  663.      LOCATE 1, 16: PRINT "Decimal";
  664.      LOCATE 2, 15: PRINT LEFT$(STR$(CPY) + " ", 8);
  665.      LOCATE 3, 15: PRINT LEFT$(STR$(CPX) + "   ", 8);
  666.      LINE (0, 0)-(178, 42 * Yfactr), 12, B'Box around it
  667.      LINE (0, 0)-(116, 42 * Yfactr), 12, B'Box around it
  668.      LET a$ = "": LET B$ = "": RETURN
  669.  
  670. HELP: CLS : COLOR 15, 1: LINE (0, 0)-(639, 18 * Yfactr), 14, BF
  671.       LOCATE 1, 20: PRINT " MAPFIX.bas HELP SCREEN Ver "; Ver$
  672.       LOCATE 3, 1
  673.       
  674.       PRINT " The cursor is shown in LAT/LON, map offset and decimal degrees.   The ORIGIN,"
  675.       PRINT " CENTER and BORDER are shown (but only the CENTER and RANGE in MAPLIST.map are"
  676.       PRINT " actually used by APRS.  Labels are right justified to the point just after the"
  677.       PRINT " last letter.  CALLS & OBJECT names will be left justified."
  678.       PRINT ""
  679.       PRINT " OPERATIONS          MAP FUNCTIONS         @N - NEW FEATURE   LABEL COMMANDS"
  680.       PRINT " H - HELP SCREENS    @C- Change CENTER     @A - ADD point     @S - SHOW labels"
  681.       PRINT "^S - SAVE MAP!!!     D - map DIRECTORY     @D - DELETE point  @L - add a LABEL"
  682.       PRINT " R - RESET pointers  M - MAPLIST.apr file  @K - Kill feature   L - LABELS off"
  683.       PRINT " Q - QUIT            O - OTHER map bordrs  @M - MOVE point"
  684.       PRINT "                     @R- set map RANGE     @T - TRIM borders"
  685.       PRINT "                                               "
  686.       PRINT " MAP COMMANDS        POINTER MOVEMENTS   USGS CD ROM CMDS   DIGITIZER & GPS "
  687.       PRINT " SPACE to draw map   N - Next Feature    B - BOX PPD area   @O- OPEN COMMS"
  688.       PRINT " ARROWS  (shft)      P - Prev Feature    U - USGS overlay   @B- BEGIN new line"
  689.       PRINT " PgUP/DN (ctrl)      G - Go to Pointer   T - Test smoother "
  690.       PRINT " HOME to Cursor      F - Find point      @S- SMOOTH file"
  691.       PRINT " END to map center   +/- move POINTER    @U- USGS BUILD!    @G- GPS OVERLAY"
  692.       PRINT "                                         @J- JOIN lines"
  693.       PRINT "                                         ^R- REDRAW on/off"
  694.       PRINT : LINE (0, 190 * Yfactr)-(639, 190 * Yfactr), 15
  695.               
  696.       IF Display$ <> "HELP" THEN
  697.        LOCATE 25, 1
  698.        PRINT " HIT H AGAIN FOR MORE HELP SCREENS, OR SPACE BAR FOR MAP...";
  699.       END IF
  700.       Display$ = "HELP"
  701.       LINE (0, 0)-(634, 348 * Yfactr), 15, B
  702.       RETURN
  703.  
  704.  
  705.  
  706. REM ************* HERE IS THE CODE BROUGHT IN FROM APRS  ***************
  707.  
  708. LdMapLst: GOSUB BoxLine23: INPUT "FileSpec for MAPLIST.apr if not \APRS\MAPLIST.APR"; a$
  709.           IF a$ <> "" THEN F$ = a$ ELSE F$ = "\aprs\Maplist.apr"
  710.           OPEN F$ FOR INPUT AS #3: IF Fault <> 0 THEN RETURN
  711.           i = 1: NumGood = 0
  712.           INPUT #3, DfltY: LINE INPUT #3, a$
  713.           INPUT #3, DfltX: LINE INPUT #3, a$
  714.           INPUT #3, BestRng: LINE INPUT #3, a$: DfltR = BestRng
  715.           INPUT #3, GMToffset: LINE INPUT #3, a$
  716.     WHILE a$ <> "* BEGIN *": LINE INPUT #3, a$: WEND ' Skip comment block
  717.           REM RS = BestRng: REM center display
  718.           REM RS = 2 ^ INT(LOG(RS) / LOG(2))'Rng is intgr power of 2
  719.           REM CPX = CDX: CPY = CDY 'Cursor Posn to Center of Display
  720.     WHILE NOT EOF(3) AND i <= UBOUND(MapName$)
  721.           INPUT #3, MapName$(i), LATcen(i), LONcen(i), MapMax(i)
  722.           LINE INPUT #3, Comment$(i)' IGNORE ALL comment fields
  723.           REM now ignore maps that start with a *
  724.           IF LEFT$(MapName$(i), 1) <> "*" THEN NumGood = NumGood + 1
  725.           NumMaps = i: i = i + 1
  726.     WEND: CLOSE #3
  727.     IF NumGood >= MaxNumMAPS - 1 THEN
  728.           CLS : LOCATE 2, 5
  729.           PRINT "WARNING: Too many ACTIVE MAPS (more than"; MaxNumMAPS; ") in MAPLIST.map file for APRS"
  730.           LOCATE 4, 10: PRINT "Use EDITOR to suppress mapnames with an (*) that you don't need."
  731.           PRINT : PRINT : PRINT : MapListLoaded = -1
  732.           INPUT "HIT RETURN to continue"; a$
  733.     END IF
  734.     RETURN
  735.  
  736. ListMAPlist: IF NOT MapListLoaded THEN GOSUB LdMapLst
  737.    GOSUB ListHeader
  738.    FOR i = 1 TO NumMaps
  739.        IF i / 19 = INT(i / 19) THEN
  740.           LOCATE 25, 1: PRINT "HIT RETURN to continue"; : INPUT a$
  741.           GOSUB ListHeader
  742.        END IF
  743.        PRINT MapName$(i); TAB(14);
  744.        PRINT INT(LATcen(i) * 100) / 100; TAB(21); INT(LONcen(i) * 100) / 100;
  745.        PRINT TAB(29); MapMax(i); TAB(36); LEFT$(LTRIM$(Comment$(i)), 43)
  746.    NEXT i
  747.    
  748.    LOCATE 25, 1: PRINT "LIST COMPLETE. CONTINUE WITH NEXT MAPFIX COMMAND...";
  749.    RETURN
  750.  
  751. ListHeader: CLS
  752.    PRINT "MAPS in MAPLIST.map (*MAPS are suppressed)     [For now, use EDITOR to modify]"
  753.    PRINT :
  754.    PRINT "MAP NAME      LATcen LONcen  RANGE COmments"
  755.    PRINT "------------  ------ ------- ----- -------------------------------------------"
  756.    RETURN
  757.  
  758. DrwAndShow: IF NOT MapListLoaded THEN GOSUB LdMapLst
  759.             Display$ = "SHOW": GOSUB DrwMPaCur
  760.  
  761. ShowMaps: IF MapSize > RS / 2 THEN MapSize = RS / 2
  762.     LOCATE 25, 1: PRINT " Drawing all maps >"; MapSize;
  763.     PRINT "mi.  F3 to see smaller, F4 for bigger, SPACE to cancel.";
  764.     LINE (0, 336 * Yfactr)-(639, 349 * Yfactr), 14, B
  765.     FOR i = 1 TO NumMaps
  766.    
  767.     x = 320 + Sfac * (CDX - LONcen(i)) * Hfac
  768.     y = Ycen + Sfac * (CDY - LATcen(i)) * Yfactr
  769.       dy = MapMax(i) * Sfac * Yfactr / 60
  770.       dx = dy * 640 / (400 * Yfactr) * Lfac
  771.       C = 15
  772.       IF MapMax(i) > 32 THEN C = 14
  773.       IF MapMax(i) > 64 THEN C = 12
  774.       IF MapMax(i) > 128 THEN C = 11
  775.       IF MapMax(i) > 256 THEN C = 13
  776.           
  777.     IF MapMax(i) > MapSize THEN
  778.        LINE (x - dx, y - dy)-(x + dx, y + dy), C, B
  779.        IF y + dy > 14 * Yfactr AND y + dy < 350 * Yfactr THEN
  780.           IF x + dx > 8 * (LEN(MapName$(i)) + 1) AND x + dx < 632 THEN
  781.              LOCATE (y + dy) / (14 * Yfactr), (x + dx) / 8 - LEN(MapName$(i))
  782.              IF MapMax(i) > RS / 4 THEN PRINT MapName$(i);
  783.           END IF
  784.        END IF
  785.     END IF
  786.     NEXT i: RETURN
  787.         
  788. REM ******* here is the code added by W7KKE for overlyaying track histoiries
  789. 'This module retrieves GPS history files so you can check the accuracy of
  790. 'the map
  791.  
  792. Hstdir: CLS : PRINT "HISTORY FILES DIRECTORY": PRINT
  793.          PRINT "To display HST files, please enter the path to your xxxxxxx.HST files."
  794.          PRINT "For example, the default '\APRS\*.HST' will show all maps in the APRS"
  795.          PRINT "directory.  Similarly '*.hst' will search your present QBasic directory."
  796.          PRINT "For any other path, enter the full file specification.": PRINT
  797.          
  798.          PRINT "Enter Filespec for searching the DIRECTORY (\aprs\*.hst)";
  799.          INPUT F$: IF F$ = "" THEN F$ = "\aprs\*.hst"
  800.          IF INSTR(F$, ".") = 0 THEN F$ = F$ + ".HST"
  801.          PRINT : PRINT : FILES F$
  802.          RETURN
  803.  
  804. LoadHst: GOSUB BoxLine23
  805.        INPUT "Which history file to load (ENTER for list, Q to quit)"; F$
  806.        IF UCASE$(F$) = "Q" THEN RETURN
  807.        IF F$ = "" THEN GOSUB Hstdir: GOTO LoadHst
  808.        a = INSTR(3, F$, "."): IF a = 0 THEN F$ = F$ + ".hst"
  809.        Fault = 0: F$ = UCASE$(F$): OPEN F$ FOR INPUT AS #3
  810.    
  811.        IF Fault = 53 OR Fault = 62 THEN Fault = 0: RETURN
  812.        GOSUB BoxLine23: PRINT "Loading track history from "; F$
  813.        
  814.     DO WHILE NOT EOF(3)
  815.        i = i + 1
  816.        INPUT #3, a$
  817.        HLAT(i) = VAL(MID$(a$, 26, 2)) + (VAL(MID$(a$, 28, 5)) / 60)
  818.        HLONG(i) = VAL(MID$(a$, 35, 3)) + (VAL(MID$(a$, 38, 5)) / 60)
  819.        maxhist = i
  820.      LOOP
  821.      CLOSE #3: Histloaded = -1
  822.      GOSUB BoxLine23: PRINT "File loading is complete.  GPS data is plotted."
  823.      REM fall through...
  824.  
  825. DrawHist:  'put history track on map
  826.      IF Histloaded THEN
  827.         size = 3: IF RS < 2 THEN size = size * 2 / RS
  828.         FOR i = 1 TO maxhist
  829.             HMX = 320 + KP * (CDX - HLONG(i)) * ppdV * Hfac'new
  830.             HMY = Ycen + KP * (CDY - HLAT(i)) * ppdV * Yfactr
  831.             CIRCLE (HMX, HMY), size, 13
  832.         NEXT i
  833.      END IF
  834.      RETURN
  835.  
  836. NewMap: CLS : PRINT "BEGINNING A NEW MAP FROM SCRATCH...": PRINT
  837.    PRINT "All points in an APRS map are measured as an"
  838.    PRINT "offset to the right and down from an origin."
  839.    PRINT
  840.    INPUT "Enter the LATITUDE  of the ORIGIN in decimal degrees"; LATo
  841.    INPUT "Enter the LONGITUDE of the ORIGIN in decimal degrees"; LONo
  842.    PRINT
  843.    PRINT "Choose the number of pixels per degree to set the map scale:"
  844.    PRINT
  845.    PRINT "Approximate size    Range from center  Pixels/Deg"
  846.    PRINT "----------------    -----------------  ----------"
  847.    PRINT "Big state or region        250             120"
  848.    PRINT "Typical state              100             300"
  849.    PRINT "Several Counties            50             600"
  850.    PRINT "Typical VHF range           25            1200"
  851.    PRINT "City streets (7.5 min maps) 12            2400"
  852.    PRINT
  853.    INPUT "Enter desired Pixels/Deg"; ppdB
  854.    IF ppdB = 0 THEN GOTO NewMap
  855.    REM In following lines, 500 is half of 999 (maximum nominal value for pts)
  856.    LATcen = LATo - (500 * Yfactr / ppdB)
  857.    LONcen = LONo - (500 / ppdB)
  858.    GOSUB StartMap: ppdV = ppdB
  859.    CLS : PRINT "YOU ARE NOW READY TO DRAW A NEW MAP...": PRINT : PRINT
  860.    PRINT "A white border has been drawn around the maximum size permitted for this map"
  861.    PRINT
  862.    PRINT "USING CURSOR WITHOUT DIGITIZER:  Move coursor to starting point for a NEW"
  863.    PRINT "feature and hit ALT-N.  Then enter new feature name (for reference purposes)"
  864.    PRINT "and continue moving cursor to the next point and hit ALT-A to add more points."
  865.    PRINT "Continue in this fashion, using ALT-N whenever you want to begin a NEW feature."
  866.    PRINT
  867.    PRINT "USING A DIGITIZER:  First, use ALT-O once to OPEN the digitizer COM port.  Then"
  868.    PRINT "use ALT-B to BEGIN each new map feature.  Enter the name and color of the new"
  869.    PRINT "feature.  Then use the digitizer mouse to add more points."
  870.    PRINT : PRINT : PRINT
  871.    PRINT "Add LABELS on the map at the current cursor location by using the ALT-L key. "
  872.    PRINT
  873.    PRINT "When you are finished, be sure to SAVE the map using the CTRL-S command..."
  874.    PRINT : PRINT : PRINT
  875.    PRINT "FOR HELP, REMEMBER THE  H  KEY!"
  876.    PRINT : PRINT : PRINT "Hit ENTER to proceed..."; : INPUT a$
  877.    RETURN
  878.  
  879. StartMap: REM This called by NEW and in middle of USGS build
  880.    LATcen$ = "LAT of CENTER": LONcen$ = "LON of CENTER"
  881.    MapRng = 60 * 500 * Yfactr / ppdB: REM 500 is half of full map size
  882.    MapRng$ = "Map range from center"
  883.    VS$ = "Pixels per degree"
  884.    MinRng = 1: MR$ = "No longer used"
  885.    TextLine$ = "NEW Map generated by MAPFIX.bas routine..."
  886.    IF Key$ = "NEW" THEN RS = 2 ^ INT(LOG(MapRng) / LOG(2))'Rng is intgr of VERTrng
  887.    CDX = LONcen: CDY = LATcen: CPX = CDX: CPY = CDY
  888.    nmp = 1: nml = 0
  889.    LNi = 1: LN$(1) = "Labels begin here"
  890.    RETURN
  891.  
  892.  
  893. DigiInit: CLS : PRINT : Digitizer = -1
  894.    PRINT "This routine will replace many CURSOR functions with the Digitizer's MOUSE."
  895.    PRINT "Assuming your digitizer can output an X,Y,C format."
  896.    PRINT
  897.    PRINT "Only Mercator projection charts will give absolutely accurate results.  Other"
  898.    PRINT "types, Lambert Conformal, Conical, etc will induce distortions."
  899.    PRINT
  900.    PRINT "It has not been tested with East Longitude or South Latitude."
  901.    PRINT : PRINT
  902.    PRINT "The digitizr should operate at 9600,N,8,1 in POINT mode with 200 LPI resolution."
  903.    PRINT "The FORMAT outputs X,Y,C values separated by commas (C is for button pressed."
  904.    PRINT
  905.    PRINT "Set up the digitizer according to your model's instructions.  For the model"
  906.    PRINT "23360, use the drawing board menu by pressing the mouse button 0 on the SETUP"
  907.    PRINT "label so that the LED is ON.  Then move the mouse to each other label and"
  908.    PRINT "use the 0 button to toggle the value ON or off as follows:"
  909.    PRINT
  910.    PRINT "POINT is ON                             PARITY 7/8 and 1 are ON "
  911.    PRINT "BAUDRATE 3 is ON                        FORMAT is ON ON off ON"
  912.    PRINT "DATA RATE doesn't matter                RESOLUTION off off ON"
  913.    PRINT : PRINT
  914.    INPUT "Is DIGITIZER connected to COM1 or COM2 (1)"; a$
  915.    IF a$ <> "2" THEN a$ = "COM1" ELSE a$ = "COM2"
  916.  
  917.    Port$ = a$ + ":9600,N,8,1,cs0,ds0,cd0"
  918.    OPEN Port$ FOR RANDOM AS #1
  919.    
  920.  
  921.    CLS : PRINT "FIRST LETS TEST THE DIGITIZER, AND GET THE MAP ON STRAIGHT.": PRINT
  922.    PRINT "Move your mouse (or pen) and hit the 0 button (or touch tablet) to see if the"
  923.    PRINT "digitizer is outputting in the desired format.  While doing this, it is a good"
  924.    PRINT "idea to verify that your map is on straight.  The Y values from the mouse"
  925.    PRINT "should give the same values for the same LATITUDE line on both the right and"
  926.    PRINT "left edges of the map.  If not, move your map to get it horizontal."
  927.    PRINT
  928.    PRINT "OUTPUT FORMAT:"
  929.    PRINT
  930.    PRINT "XXXXX,YYYYY,APn (Only the X and Y values are used (4 or 5 digits is ok)"
  931.    PRINT
  932.    LOCATE 25, 1: PRINT "Hit ENTER and press 0 button on mouse to continue...";
  933.    LOCATE 13, 1
  934.    DO UNTIL INKEY$ <> "": LINE INPUT #1, a$: PRINT a$: LOOP
  935.   
  936.    CLS : PRINT
  937.    PRINT "NEXT YOU MUST ESTABLISH THE SCALE OF YOUR DIGITIZER."
  938.    PRINT
  939.    PRINT "The scale is established by two points, the first near the"
  940.    PRINT "upper left corner, the second near the lower right corner."
  941.    PRINT
  942.    PRINT "To get the best accuracy on maps not exactly MERCATOR, use points within the "
  943.    PRINT "area where you are working, not on the extreme corners.  IE:  choose points"
  944.    PRINT "that are in the center of the upper left quadrant and the lower right quadrant."
  945.    PRINT
  946.    PRINT "To establish the upper left reference point:"
  947.    INPUT "     Enter lat  (deg,min)"; LATref1, M: LATref1 = LATref1 + M / 60
  948.    INPUT "     Enter long (deg,min)"; LONref1, M: LONref1 = LONref1 + M / 60
  949.    PRINT
  950.  
  951.    PRINT "Place the mouse on the upper left point and press the 0 button."
  952.    LINE INPUT #1, a$: SOUND 150, 3
  953.         digix1 = 5000 - VAL(LEFT$(a$, 5))
  954.         digiy1 = VAL(MID$(a$, 7, 5))
  955.         PRINT "Digitizer reads "; digix1, digiy1; " for this point.": PRINT
  956.  
  957.    PRINT "NOW Establish the lower right reference point:"
  958.    INPUT "      Enter lat  (deg,min)"; LATref2, M: LATref2 = LATref2 + M / 60
  959.    INPUT "      Enter long (deg,min)"; LONref2, M: LONref2 = LONref2 + M / 60
  960.    PRINT
  961.    PRINT "Place digitizer pen on lower right point."
  962.    LINE INPUT #1, a$: SOUND 150, 3
  963.         digix2 = 5000 - VAL(LEFT$(a$, 5))
  964.         digiy2 = VAL(MID$(a$, 7, 5))
  965.         PRINT "Digitizer reads "; digix2, digiy2; " for this point.": PRINT
  966.  
  967.    REM Find delta lat/long between reference points
  968.    REM Calculate degrees per x/y unit
  969.         degx# = (LONref1 - LONref2) / (digix1 - digix2)
  970.         degy# = (LATref1 - LATref2) / (digiy1 - digiy2)
  971.  
  972.    CLS : PRINT "YOU ARE NOW READY TO USE THE DIGITIZER TO ENTER POINTS INTO MAPFIX..."
  973.    PRINT
  974.    PRINT "The digitizer works just about like the cursor and arrow keys in MAPFIX.  Any"
  975.    PRINT "point identified by the digitizer will be ADDED just as if you had hit ALT-A."
  976.    PRINT "All points are added to a feature after the current MapPoint identified by the"
  977.    PRINT "white circle. "
  978.    PRINT
  979.    PRINT "With the digitizer, do NOT use the ALT-N NEW command which always begins at the"
  980.    PRINT "current cursor location.  For the digitizer, use ALT-B to BEGIN a new feature."
  981.    PRINT "You will be asked to identify the name and color of the new feature.  From then"
  982.    PRINT "on, just move the digitizer mouse (or pen) to ADD new points.  "
  983.    PRINT
  984.    PRINT "If your digitizer mouse has 4 buttons, use the first (left) button for ADDing "
  985.    PRINT "points, use the 4th (right) button to just move the cursor with no action."
  986.    PRINT : PRINT
  987.    INPUT "Hit ENTER to continue with MAPFIX..."; a$
  988.    RETURN
  989.  
  990. GetXY: LINE INPUT #1, a$: SOUND 150, 3
  991.        a = INSTR(a$, ","): IF a = 0 THEN RETURN
  992.            x = 5000 - VAL(LEFT$(a$, a - 1))
  993.        B = INSTR(a + 1, a$, ","): IF B = 0 THEN B = LEN(a$)
  994.        y = VAL(MID$(a$, a + 1, B - (a)))
  995.        Btn = VAL(RIGHT$(a$, 1))
  996.        CPY = ((y - digiy2) * degy#) + LATref2
  997.        CPX = ((x - digix2) * degx#) + LONref2
  998.        IF LOC(1) <> 0 THEN a$ = INPUT$(LOC(1), #1)'Clear input buffer
  999.        RETURN
  1000.  
  1001. BoxPPD: GOSUB BoxLine23: INPUT "Enter the desired PPD"; a$: ppdB = VAL(a$)
  1002.    dy = (30000 / ppdB) * Sfac * Yfactr / 60
  1003.    dx = dy * 640 / (400 * Yfactr) * Lfac
  1004.    GOSUB BoxLine23: LINE (CUX - dx, CUY - dy)-(CUX + dx, CUY + dy), 13, B
  1005.    PRINT "The box represents the largest APRS map that can be made with that scale."
  1006.    RETURN
  1007.         
  1008.  
  1009. GetUSGS: REM This used for both U=OVERLAY and by ALT-U = USGS BUILD!
  1010.    USGS = -1: ni = 0: nt = 0: j = 0: NumLines = 0
  1011.    IF Key$ <> "U" THEN
  1012.       IF ppdB <> 0 THEN ppdV = ppdB
  1013.       LATo = CDY + (500 * Yfactr / ppdV)
  1014.       LONo = CDX + (500 / ppdV)
  1015.       GOSUB BoxLine23: PRINT "Improve LAT ORIGIN of "; LATo; : INPUT LATo
  1016.       GOSUB BoxLine23: PRINT "Improve LON ORIGIN of "; LONo; : INPUT LONo
  1017.       dx = ppdV * (LONo - CDX)
  1018.       dy = ppdV * (LATo - CDY)
  1019.       KP = 100 * 100 / (RS * ppdV)
  1020.       LATcen = CDY: LONcen = CDX: GOSUB StartMap
  1021.       LATtext$ = "Decimal LATITUDE  of map ORIGIN"
  1022.       LONtext$ = "Decimal LONGITUDE of map ORIGIN"
  1023.       TextLine$ = "Map generated by MAPFIX from USGS 2,000,000:1 CD ROM (data valid mid-1980's)"
  1024.       END IF
  1025.    GOSUB BoxLine23: INPUT "Enter path and filename of XTRACTED USGS data file"; a$
  1026.    a = INSTR(a$, "."): IF a = 0 THEN a$ = a$ + ".GRF"
  1027.    OPEN a$ FOR INPUT AS #3
  1028.    IF Fault <> 0 THEN RETURN
  1029.    GOSUB BoxLine23: INPUT "Enter RANK cutoff. (all roads = 99) for ST and WB use 20, or 10"; a$
  1030.    MaxRnk = VAL(a$)
  1031.    REM PRINT "raw data format.....", "   LineID", "#-Rnk-Atbts", "  NumPts"
  1032.    LOCATE 5, 67: PRINT "RANGE: "; INT(30000 / ppdV): LOCATE 24, 1
  1033.    IF Key$ = "U" THEN
  1034.       PRINT "While USGS OVERLAYED, do not redraw map or you will have to do it again...";
  1035.    ELSE PRINT "Blue circles start lines, Green Box ends.  Red points discarded, Yellow Kept!";
  1036.    END IF
  1037.  DO UNTIL EOF(3) OR LNi = MaxNumLines - 1
  1038.    NumLines = NumLines + 1
  1039.    LOCATE 1, 61: PRINT "TOTAL POINTS: "; nt
  1040.    LOCATE 2, 61: PRINT "POINTS USED:  "; ni
  1041.    LOCATE 3, 61: PRINT "TOTAL LINEs:  "; NumLines
  1042.    LOCATE 4, 61: PRINT "LINEs USED:   "; LNi
  1043.    a$ = INPUT$(20, 3): REM PRINT a$;
  1044.         LnID$ = LEFT$(a$, 7)
  1045.         Rank$ = MID$(a$, 8, 2): Rank = VAL(Rank$): LOCATE 6, 67: PRINT "RANK: "; Rank
  1046.         Npts$ = MID$(a$, 10, 6): Npts = VAL(Npts$)
  1047.         AtCd$ = MID$(a$, 16, 5)
  1048.         a$ = LTRIM$(LnID$) + "-" + Rank$ + "-" + AtCd$
  1049.         REM PRINT , LnID$, a$, Npts$
  1050.    IF Rank < 24 THEN SavClr = 4 ELSE SavClr = 7
  1051.    IF Rank < 20 THEN SavClr = 12
  1052.    IF Rank < 14 THEN SavClr = 10
  1053.    LineOK = 0: IF Key$ <> "U" THEN GOSUB BeginF
  1054.    FOR i = 1 TO Npts
  1055.     a$ = INPUT$(20, 3): IF Rank > MaxRnk THEN GOTO Skp
  1056.     REM IF VAL(Rank$) > 99 THEN GOTO Skp
  1057.     LA = VAL(LEFT$(a$, 2)) + VAL(MID$(a$, 3, 2)) / 60 + VAL(MID$(a$, 5, 2)) / 3600
  1058.     LO = VAL(MID$(a$, 8, 3)) + VAL(MID$(a$, 11, 2)) / 60 + VAL(MID$(a$, 13, 2)) / 3600
  1059.     IF Key$ = "U" THEN
  1060.         REM Following lines used to limit points if just doing an OVERLAY only
  1061.         IF LA > CDY + RS / 60 OR LA < CDY - RS / 50 THEN GOTO Skp 'off screen
  1062.         IF LO > CDX + RS / 40 OR LO < CDX - RS / 40 THEN GOTO Skp
  1063.         REM s$ = MID$(a$, 16, 5)
  1064.         REM PRINT S$, LA, LO
  1065.         END IF
  1066.     y% = (LATo - LA) * ppdV: x% = (LONo - LO) * ppdV
  1067.     IF x% > 1100 OR x% < -100 THEN GOTO Skp 'this ignores points off PPD
  1068.     IF y% > 1100 OR y% < -100 THEN GOTO Skp 'scale
  1069.     LineOK = -1
  1070.        X1 = 320 + KP * (x% - dx) * Hfac
  1071.        Y1 = Ycen + KP * (y% - dy) * Yfactr
  1072.        IF i > 2 THEN
  1073.           REM LINE (x, y)-(X1, Y1), 6
  1074.           dd = LO - LOb: IF dd = 0 THEN dd = .0000001
  1075.           dn = LA - LAb
  1076.           s = dn / dd' Note that 1>s>.01 for Xdelta of 1 to 100
  1077.           IF ABS(s) < .1 AND ABS(Lsp) < .1 THEN
  1078.              sd = 1
  1079.           ELSEIF ABS(s) > 10 AND ABS(Lsp) > 10 THEN sd = 1
  1080.           ELSEIF ABS(dd) < .004 AND ABS(dn) < .004 THEN sd = 1
  1081.           ELSEIF s <> 0 THEN sd = Lsp / s
  1082.           ELSE sd = 0
  1083.           END IF
  1084.           IF sd > 2 OR sd < .5 OR i = Npts THEN
  1085.              CIRCLE (x, y), 1, 14
  1086.              IF Key$ <> "U" THEN GOSUB KeepLine
  1087.           ELSE CIRCLE (x, y), 1, 4
  1088.           END IF
  1089.           Lsp = s: nt = nt + 1
  1090.        ELSE Lsp = 0: CIRCLE (X1, Y1), 4, 9
  1091.           IF Key$ <> "U" THEN GOSUB KeepLine 'keeps first two lines
  1092.        END IF
  1093.        LAb = LA: LOb = LO
  1094.        x = X1: y = Y1
  1095.        
  1096. Skp: NEXT i
  1097.      IF Key$ <> "U" THEN
  1098.         IF LineOK THEN nmp = nmp - 1: Z = Z - 1: ni = ni + 1 ELSE GOSUB CanclF
  1099.      END IF
  1100.      LINE (x - 3, y - 3)-(x + 3, y + 3), 10, B ' Last Point
  1101.  LOOP
  1102.  IF LNi > MaxNumLines - 2 THEN LOCATE 12, 20: PRINT "PROCESSING STOPPED... TOO MANY LLINES!..."
  1103.    CLOSE #3
  1104.    RETURN
  1105.  
  1106. KeepLine: x%(Z) = x%: y%(Z) = y%: nmp = nmp + 1: Z = Z + 1: ni = ni + 1: RETURN
  1107.  
  1108.  
  1109. Scrunch: i = 0: Pt = 0: nt = 0: ni = 0: GOSUB BoxLine23
  1110.     INPUT "Enter slope filter ratio 1.2 to 5 (typically 1.5)"; a$
  1111.     IF a$ = "" THEN slope = 1.5 ELSE slope = VAL(a$)
  1112.     DO UNTIL i >= nmp - 1
  1113.        i = i + 1
  1114.        X1 = 320 + KP * (x%(i) - dx) * Hfac
  1115.        Y1 = Ycen + KP * (y%(i) - dy) * Yfactr
  1116.        IF x%(i) <> 0 THEN
  1117.           Pt = Pt + 1
  1118.           IF Pt > 2 THEN
  1119.              LINE (x, y)-(X1, Y1), 6
  1120.              dd = x - X1: IF dd = 0 THEN dd = .01
  1121.              dn = y - Y1
  1122.              dst = ((dd * dd) + (dn * dn)) ^ .5
  1123.              s = dn / dd' Note that 1>s>.01 for Xdelta of 1 to 100
  1124.              IF s = 0 THEN s = .05
  1125.              IF ABS(s) < .2 THEN s = .2 * SGN(s)
  1126.              IF ABS(s) > 5 THEN s = 5 * SGN(s)
  1127.              IF ABS(s) <= .2 AND ABS(Lsp) <= .2 THEN
  1128.                 sd = 1
  1129.              ELSEIF ABS(s) >= 5 AND ABS(Lsp) >= 5 THEN sd = 1
  1130.              ELSE sd = Lsp / s
  1131.              END IF
  1132.              IF ABS(dd) > 50 * KP OR ABS(dn) > 30 * KP THEN sd = 0
  1133.              REM IF ABS(dd) < 5 OR ABS(dn) < 4 THEN sd = 1
  1134.              IF sd > slope OR sd < 1 / slope OR x%(i + 1) = 0 OR NumRej > 4 THEN
  1135.                 ni = ni + 1: CIRCLE (x, y), 2, 14: NumRej = 0
  1136.                 
  1137.              ELSE CIRCLE (x, y), 1, 4: NumRej = NumRej + 1
  1138.                 IF Key$ <> "T" THEN
  1139.                    i = i - 1: nmp = nmp - 1
  1140.                    FOR ii = i TO nmp
  1141.                        x%(ii) = x%(ii + 1): y%(ii) = y%(ii + 1)
  1142.                    NEXT ii
  1143.                 END IF
  1144.              END IF
  1145.              Lsp = s: nt = nt + 1
  1146.           ELSE Lsp = 0: nt = nt + 1: ni = ni + 1: CIRCLE (X1, Y1), 4, 9
  1147.           END IF
  1148.        ELSE Pt = 0: nt = nt + 1: ni = ni + 1
  1149.           LOCATE 1, 61: PRINT "TOTAL POINTS: "; nt
  1150.           LOCATE 2, 61: PRINT "SAVED POINTS: "; ni
  1151.        END IF
  1152.        x = X1: y = Y1
  1153.     LOOP
  1154.     RETURN
  1155.  
  1156. Join: REM Search for end=begin point values and CONCATONATE if equal!
  1157.    LNptr = 0: i = 0: k = 0: GOSUB BoxLine23: PRINT "Lines joined: ";
  1158.    DO UNTIL i >= nmp
  1159.       i = i + 1
  1160.        IF x%(i) = x%(i + 2) AND y%(i) = y%(i + 2) AND y%(i + 1) = LColor THEN
  1161.           nmp = nmp - 2: LNi = LNi - 1: k = k + 1: LOCATE 23, 15: PRINT k
  1162.           FOR j = i + 1 TO nmp: x%(j) = x%(j + 2): y%(j) = y%(j + 2): NEXT j
  1163.           FOR j = LNptr TO LNi: LN$(j) = LN$(j + 1): NEXT j
  1164.        ELSEIF x%(i) = 0 THEN LColor = y%(i): LNptr = LNptr + 1
  1165.        END IF
  1166.    LOOP: GOTO DrawMap
  1167.  
  1168. END
  1169.  
  1170.